AEU Int. J. Electron. Commun. 
58 (2004) No. 1, 1-3 



© Gustav Fischer Verlag 
Jena 



1 



Size-Dependent Bruggeman Approach for Dielectric-Magnetic 

Composite Materials 

Akhlesh Lakhtakia, Tom G. Mackay 



Abstract Expressions arising from the Bruggeman approach 
for the homogenization of dielectric-magnetic composite mate- 
rials, without ignoring the sizes of the spherical particles, are 
presented. These expressions exhibit the proper limit behavior. 
The incorporation of size dependence is directly responsible for 
the emergence of dielectric-magnetic coupling in the estimated 
relative permittivity and permeability of the homogenized com- 
posite material. 
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1. Introduction 

The objective of this communication is to introduce a 
size-dependent variant of the celebrated Bruggeman ap- 
proach [l] Eq. 32] and thereby couple the dielectric and 
magnetic properties of a particulate composite material 
(PCM) with isotropic dielectric-magnetic constituent ma- 
terials. 

Homogenization of PCMs has been a continuing theme 
in electromagnetism for about two centuries 1 2 1 . The most 
popular approaches consider theparticles to be vanish- 
ingly small, point-like entities J3]E]. Much of the liter- 
ature is devoted to dielectric PCMs |3 5|, with applica- 
tion to magnetic PCMs following as a result of electro- 
magnetic duality Sec. 4-2.3]. When PCMs with both 
dielectric and magnetic properties are considered, no cou- 
pling arises between the two types of constitutive proper- 
ties if the particles are vanishingly small. It is this coupling 
that has gained importance in the last few years, with the 
emergence of metamaterials Q. 

Investigation of scattering literature quickly reveals 
that dielectric-magnetic coupling in PCMs emerges only 
when particles are explicitly taken to be of nonzero size 
E |9| 1101 . although the particle size must still be elec- 
trically small for the concept of homogenization to re- 
main valid p. xiii], II II . To the best of our knowl- 
edge, available homogenization formulas for dielectric- 
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magnetic PCMs that also account for dielectric-magnetic 
coupling are applicable only to dilute composites because 
they are set up using the Mossotti-Clausius approach (also 
called the Lorenz-Lorentz approach and the Maxwell Gar- 
nett approach 11 21 ). Use of the Bruggeman approach is 
preferred, while maintaining the particle size as nonzero, 
for nondilute composites 1 12 1. 

Accordingly, in Section [2] we apply the Bruggeman 
approach to derive size-dependent homogenization for- 
mulas for dielectric-magnetic PCMs comprising spher- 
ical particles. Sample results are discussed and conclu- 
sions are drawn therefrom in Section [3] An cxp(— iu>t) 
time-dependence is implicit, with u> being the angular fre- 
quency. The free-space (i.e., vacuum) wavenumber is de- 
noted by ko. 



2. Theory 

Let us consider a particulate composite material with L 
constituent materials. The relative permittivity and the rel- 
ative permeability of the £th constituent material, I £ 
[1, L], are denoted respectively by q and tii, the radius 
of the spherical particles of that material is denoted by Rg, 
and the volumetric fraction by fg. Clearly, 



(i) 



Our task is to estimate chcm and hhcm, which are the 
relative permittivity and the relative permeability of the 
homogenized composite material (HCM). 

According to the Bruggeman approach |4 10 1, the fol- 
lowing two equations have to be solved: 



'Br 



(2) 



£=1 



Here, a 



a/b 



and 



a 



a/b 



are the polarizability density and 



the magnetizability density, respectively, of an electrically 
small sphere of material a embedded in material b. In the 
limit of the particulate radius tending to zero, expressions 
of these two densities are available as follows [ 13 1: 



a/b n 



£a+2et, 
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3(l b 



-Mb 



(3) 



However, when the sphere radius is nonzero, the foregoing 
expressions mutate to include both the radius R a of the 
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embedded sphere and the refractive index 



rib = V£6M6 
of the embedding material; thus 1 10 1 

a/b 



a /b _ 



ot * e*(l-2T a/b )+2e b (l+T a/b ) 
°^ AIa (l-2T £l/6 )+2 Al6 (l+r £l/[ ,) 



where 



fyft = (1 - ikoR a n b ) exp(ik R a n b ) - 1 



(4) 



(5) 



(6) 



More complicated expressions than @ can be devised by 
using the Lorenz-Mie-Debye formulation for scattering 
by a sphere |8|, but do not lead to significantly different 
results for electrically small spheres. In the limit R a — ► 0, 
expressions reduce to Q because 



lim t, 



a/b 



0. 



(7) 



Clearly, the incorporation of particle size-dependence 
via (|5jl in (0 leads to a coupling of the relative permittivi- 
ties and the relative permeabilities. 



3. Results and Conclusion 

In order to investigate the properties of 0, let us simplify 
it for a two-constituent composite material: L = 2. Ex- 
pressions © for the size-dependent Bruggeman approach 
then read as follows: 



,1 £l(l-2T 1/Bl .) + 2 £f3r (l + T 1/Br ) 



= /i 

^ e2(l-2r 2/Br )+2e Br (l+T 2/Br ) 
•' 1 e 1 (l-2r 1/Br )+2 M B r (l+T 1/Br ) 



At2(l-2r 2/Br )+2^ Br (l+r 2/Br ) , 



(8) 



These two coupled equations have to be solved together in 
order to obtain the estimates e Br and fiBr of chcm and 
/J-hcm as functions of ko, /i, i?i and R 2 . 

Equations (|8} have to be solved iteratively, and the 
Newton-Raphson method is very useful for that purpose 
II 141 Sec. 6.5.2]. Typically, this method requires an initial 
guess, which can be supplied using the Maxwell Garnett 
approach |2 3|. If /i > f 2 , then material 1 should be 
treated as the host material while material 2 is dispersed 
in particulate form; and the size-dependent Maxwell Gar- 
nett estimates of chcm and \ihcm are then obtained as 
follows: 



£mg,i = ei + (1 - fx) 



^71" 

MW0-*r 



fMG.i = m + (i - h)- 



(9) 



On the other hand, the size-dependent Maxwell Garnett 
estimates 



£JMG,2 — £2 + /l 



,1/1 



1/2 



(10) 



3^2 



appear more appropriate when f 2 > f\. Incorporation of 
size dependence couples dielectric and magnetic proper- 
ties also in and JlOi . 
Let us note that the limits 
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satisfied by the solutions of (|8} are physically correct, and 
are not affected by the incorporation of size dependence in 
the Bruggeman approach. In contrast, the size-dependent 
Maxwell Garnett expressions (|9j and (II Oi do not exhibit 
physically reasonable limits when the host material van- 
ishes; i.e., 



lim/^o 
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and analogously for £a/g,2 and (img,2- The foregoing lim- 
its are borne out by the plots of chcm and Hhcm versus 
fi presented in Figures|20 

Figure ^ presents estimates of the real and imaginary 
parts of the relative permittivity and the relative perme- 
ability of the HCM when e x = 1.5, = 1, e 2 = 5 + i0.2, 
and fi2 = 2 + iO.l, and the sizes R\ = R 2 = 0. Calcu- 
lations for the relative permittivity and the relative perme- 
ability then decouple from each other. 

The analogous plots in Figure|2]were drawn for i?i = 
i?2 7^ 0. These plots are quite different from those in 
the preceding figure. The imaginary parts of chcm and 
Hhcm appear to be more affected by the size dependence 
than the real parts. Indeed, were both constituent materials 
totally nondissipative, the imaginary parts of T a / b -terms 
would still give rise to imaginary parts of both chcm and 
Uhcm 11121 . We also conclude from comparing Figures 
]l\ and |2] that dielectric-magnetic coupling proportionally 
affects the imaginary parts of the HCM constitutive pa- 
rameters more than their real parts. 

There is no reason for the particles of both constituent 
materials to be of the same size (or have the same dis- 
tribution of size, in general). The plots in Figure [3] were 
drawn for R 2 — 3i?i. Clearly from this figure and Fig- 
ure |2] the effect of different particle sizes on dielectric- 
magnetic coupling can be substantial. 

The permeability contrast between the two constituent 
materials chosen for Figures is l ess than the permit- 
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Fig. 1. Estimates of the real and imaginary parts of the relative 
permittivity and the relative permeability of a homogenized com- 
posite material (HCM) with two constituent materials (ei = 1.5, 
/ii = 1, £2 = 5 + i0.2, and [1% — 2 + iO.l) as functions of 
the volumetric fraction /i = 1 — /a. Size-independent Maxwell 
Garnett approach with material 1 as the host material (dashed line); 
Size-independent Maxwell Garnett approach with material 2 as the 
host material (dotted line); Size-independent Bruggeman approach 
(solid line). fco^?i = koRi = 0. 




Fig. 2. Estimates of the real and imaginary parts of the relative 
permittivity and the relative permeability of a homogenized com- 
posite material (HCM) with two constituent materials (ei = 1.5, 
fXi = 1, €i — 5 + i0.2, and fj,2 = 2 + iO.l) as functions of the 
volume fraction /i = 1 — f-z. Size-dependent Maxwell Garnett 
approach with material 1 as the host material (dashed line); Size- 
dependent Maxwell Garnett approach with material 2 as the host 
material (dotted line); Size-dependent Bruggeman approach (solid 
line), fco^i = k R 2 = 0.2. 



tivity contrast. We notice that the effect of size depen- 
dence on [ahcm is less than on chgm- This implies that 
dielectric-magnetic coupling affects the more contrasting 
constitutive parameter more. 

To conclude, we have here implemented the Bruggeman 
approach for the homogenization of dielectric-magnetic 
composite materials, without ignoring the sizes of the 
spherical particles. These expressions exhibit the proper 
limit behavior. The incorporation of size dependence 
is directly responsible for the emergence of dielectric- 



Fig. 3. Same as Figure|2| except that fco Ri = 0.2 and ko R2 = 0.6. 



magnetic coupling in the estimated relative permittivity 
and permeability of the homogenized composite material. 
The size-dependent Bruggeman estimates are compared 
with the size-dependent Maxwell Garnet estimates, which 
do not necessarily evince the proper limit behavior and are 
therefore applicable to dilute composite materials. 
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